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Abstract 

The miltifractal properties and scaling behaviour of the exchange 
rate variations of the Iranian rial against the US dollar from a daily 
perspective is numerically investigated. For this purpose the multifrac- 
tal detrended fluctuation analysis (MF-DFA) is used. Through mul- 
tifractal analysis, the scaling exponents, generalized Hurst exponents, 
generalized fractal dimensions and singularity spectrum are derived. 
Moreover, contribution of two major sources of multifractality, that is, 
fat-tailed probability distributions and nonlinear temporal correlations 
are studied. 

Keywords: Multifractality, Scaling, Rial-dollar exchange rate, Financial 
markets. 

1 Introduction 

For more than two decades, there has been considerable interest in the in- 
vestigation of the scaling behaviour on fractal models. The pioneering work 
on fractals introduced the concept of fractals and showed some rela- 
tion between self-similar fractals and self-affine fractals. Self-affine fractals 
[U E3 IU IS] constitute random and complicated structure and have been 
applied to a broader range of problems, such as the Eden and ballistic 
deposition model |3 El E] , mountain heights, clouds, coast lines, and 
cracks. Among other examples of many fractal models, the self-avoiding 
random walk, random resistor, polymer bonds, turbulences, chaotic mo- 
tions can be mentioned [U QJ EH > e£c Specially, the real 
data from different financial markets show apparent multifractal properties 
US [T3 Eni dZl UHl Uni EDI EH Moreover, recently, financial analysis of 
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foreign exchanges has became one of the outstanding topics in econophysics 
|23| . Many of these researches apply multifractal analysis framework as 
a basic framework. In universal multifractal framework, the statistics of 
the data are fully described with some parameters, taking into account two 
complementary aspects of financial time series: the multiple scaling and the 
Pareto probability distributions, which is a generic feature of multifractal 
processes |2"1~] , 

It has been shown that there are two main factors leading to multifractal be- 
haviour of financial time series, nonlinear time correlations between present 
and past events and the heavy-tailed probability distributions of functions. 
Based on ref. |21| , for the stocks, the main contribution to multifractality 
comes from a broad distribution of returns while a long memory present 
in this kind of data contributes only marginally. It should be noted that 
the nature of correlations leading to the multifractal dynamics of the varia- 
tions is strongly nonlinear and, curiously, cannot be simply related to some 
well-known correlation type like a slowly decreasing volatility autocorrela- 
tion with an imposed daily pattern. For example, one even has to consider 
the nonlinear dependencies in the volatility itself in order to reveal how the 
temporal correlations contribute to multifractality in the stock market and 
foreign exchange data. 

In this paper, the rial-dollar exchange rate data is studied with the focus 
on their fractal properties. The multifractal detrended fluctuation analysis 
is applied which is a well-established method of detecting scaling behaviour 
of time series. In Section 1, theoretical backgrounds including MFDFA 
method, sources of multifractality, multifractality finger prints and strength 
of multifractality are reviewed. Data are described in Section 2. Numerical 
results are presented in Section 3 and finally, conclusions are given in Section 
4. 

2 Theoretical backgrounds 
2.1 Method 

Detrended fluctuation analysis (DFA) is a scaling analysis technique pro- 
viding a simple quantitative parameter-the scaling exponent a-to represent 
the correlation properties of a time series The advantage of DFA over 
many techniques are that it permits the detection of long-range correlations 
embedded in seemingly non-stationary time series, and also avoids the spu- 
rious detection of apparent long-range correlations that are an artifact of 
non-stationarity. Additionally, the advantages of DFA in computation of H 
over other techniques (for example, the Fourier transform) are: 

• inherent trends are avoided at all time scales; 
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• local correlations can be easily probed. 

To implement the DFA, let us suppose there is a time series, N(i)(i 
..,JV m(B ). The time series N(i) is integrated: 



y(j) = - (N)} (i) 

i=i 



where: 

1 * max 

w = w— E ( 2 ) 

1* max , 

Next iV(i) is broken up into X non-overlapping time intervals, J n , of equal 
size r where n = 0,1, ...if — 1 and K corresponds to the integer part of 
N max /r. In each box, the integrated time series is fitted by using a poly- 
nomial function, y po i(i), which is called the local trend. For order-/ DFA 
(DFA-1 if 1=1, DFA-2 if 1=2, etc.), the Z-order polynomial function should 
be applied for the fitting. The integrated time series y(i) is detrended in 
each box, and calculated the detrended fluctuation function: 

Y(i)=y(i)-y pol (i). (3) 

For a given box size s, the root mean square fluctuation is calculated: 
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1 N<max 

N — E t y «] 2 ( 4 ) 



The above computation is repeated for box sizes s (different scales) to pro- 
vide a relationship between F(s) and s. A power law relation between F(s) 
and s indicates the presence of scaling: F(s) ~ s a . The parameter a, 
called the scaling exponent or correlation exponent, represents the corre- 
lation properties of the signal: if a = 0.5, there is no correlation and the 
signal is an uncorrelated signal if ot < 0.5, the signal is anticorrelated; if 
a > 0.5, there are positive correlations in the signal. In the two latest cases, 
the signal can be well approximated by the fractional Brownian motion law 

For a further characterization of data it is meaningful to extend Eq. (15) by 
considering the more general fluctuation functions j^T] . Simply, it is achieved 
by averaging over all boxes to obtain the qth order fluctuation function 

W = l^T~ E (F 2 (s))i/r /q , (5) 

where, in general, the index variable q can take any real values except 
zero. If the analyzed signal develops fractal properties, the fluctuation func- 
tion reveals power-law scaling 

F a (s) ~ s h ^ (6) 
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for large s. The scaling exponents h(q) can be then obtained by observing the 
slope of log-log plots of F q vs.s. The family of the exponents h(q) describe 
the scaling of the qth order fluctuation function. For positive values of q, h(q) 
exponents describe the scaling behaviour of boxes with large fluctuations 
while those of negative values of q, describe scaling behaviour of boxes with 
small fluctuations |2H]- For stationary time series, the exponent h{2) is 
identical to the Hurst exponent. Thus the exponents h{q) are called as the 
generalized Hurst exponents [2E3- F° r monofractal time series which are 
characterized by a single exponent over all scales, h(q) is independent of q, 
whereas for a multifractal time series, h(q) varies with q. This dependence 
is considered to be a characteristic property of multifractal processes |28| . 
The h(q) obtained from MF-DFA is related to the Renyi exponent r(g) by 

qh{q) = r(q) + 1. (7) 

Therefore, another way to characterize a multifractal series is the singularity 
spectrum /(a) defined by [2Sj 



a = h(q) + qti(q), f(a) = q[a-h(q)]+l, (8) 

where h'(q) stands for the derivative of h(q) with respect to q. a is the Holder 
exponent or singularity strength which characterizes the singularities in a 
time series. The singularity spectrum /(a) describes the singularity content 
of the time series. Finally, it must be noted that h(q) is different from the 
generalized multifractal dimensions 

q — 1 q — 1 

that are used instead of r(q) in some papers. While h{q) is independent of 
q for a monofractal time series with compact support, D{q) depends on q in 
that case. 



2.2 Sources of multifractality 

Generally, there are two different types of sources for multifractality in time 
series: (i) due to different long-range temporal correlations for small and 
large fluctuations, and (ii) due to fat-tailed probability distributions of vari- 
ations. Both of them need a multitude of scaling exponents for small and 
large fluctuations. Two procedure is followed to find the contributions of 
two sources of multifractality and to indicate the multifractality strength: 
(i) shuffling, and (ii) phase randomization. Shuffling procedure preserves 
the distribution of the variations but destroys any temporal correlations. In 
fact, one can destroy the temporal correlations by randomly shuffling the 
corresponding time series of variations. What then remains are data with 
exactly the same fluctuation distributions but without memory. The shuffling 
procedure consists of the following steps 
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(i) Generate pairs (p, q) of random integer numbers (with p,q < N) where 
N is the total length of the time series to be shuffled. 

(ii) Swap entries p and q. 

(iii) Repeat two above steps for 20 N times. (This step ensures that order- 
ing of entries in the time series is fully shuffled.) 

In order to study the contribution of the fat-tailed variations on the multi- 
fractality, the surrogate data are used. In fact, the non-Gaussianity of the 
distributions can be weakened by creating the phase-randomized surrogates 
|29| , The Phase randomization steps are: 

(i) Take discrete Fourier transform of time series. 



(ii) Multiply the discrete Fourier transform of the data by random phases. 



(iii) Perform an inverse Fourier transform to create a phase randomized 
surrogates. 

Phase randomization preserves the amplitudes of the Fourier transform but 
randomizing the Fourier phases. This procedure eliminates nonlinearities, 
preserving only the linear properties of the original time series f30] . 



2.3 Multifractality finger prints 

One can see that in the whole q-range the generalized Hurst exponents h(q) 
can be fitted well by the formula 

w q qln2 V ' 

which corresponds to r{q) = — ln[a q +b q ] / ln2 . This formula can be obtained 
from a generalized binomial multifractal model |31j . Instead of choosing a 
and b, the Hurst exponent h(l) and the persistence exponent h(2) could be 
chosen. From knowledge of two moments, all the other momemts follow. 
Here the formula is used only to show that the infinite number of expo- 
nents h(q) can be described by only two independent parameters, a and b. 
These two parameters can then be regarded as multifractal finger prints for 
a considered time series. 



2.4 Strength of multifractality 

In the generalized binomial multifractal model, the strength of the multi- 
fractality of a time series can be characterized by the difference between the 
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maximum and minimum values of a, a max — a m i n . When q^j^- approaches 
zero for q approaching ±00, then Aa = a max — a m i n is simply given by 

Aa = h(—oo) — /i(oo) = — . (11) 

It must be noted that this parameter is identical to the width of the singu- 
larity spectrum f{ot) at / = 0. The wider singularity spectrum the richer 
multifractality. 



3 Data Analysis 

The data which is analyzed are the time series of the daily closing exchange 
rate logarithmic variations (that is, ln(P(t))/ ln(P(t+l)) for the time period 
24th September 1989, to 15th November 2003. So that our database con- 
sists of 4369 exchange rates and 4368 daily variations. The sources of this 
data is the central bank of the islamic republic of Iran. In Fig. 1 a time 
series corresponding to daily values of the Iranian rial-US dollar exchange 
rates in mentioned period is presented. A great increment in dollar price is 
seen about 16th May 1995 because of Iranian government decision on lift- 
ing the ban on foreign exchanges price variations. The Iranian government 
has managed to keep the exchange rate stable at around 8000 rials per US 
dollar ever since March 2000. Also, Table 1 provides summary statistics of 
logarithmic variations of exchange rates. According to data in Table 1, a 



Table 1: Mean, standard deviation, skewness, and kurtosis of rial-dollar 
exchange rate variations. 



Mean 


Std.Dev. 


Skewness Kurtosis 


0.00055 


0.0117 


-1.2504 49.925 



negatively large skew is seen. The probability distribution function of varia- 
tions also show a high degree of peakedness and fat tails relative to a normal 
distribution. Thus there is a clear departure from Gaussian normality. The 
departure from a Gaussian Cumulative Distribution Function (CDF) can be 
clearly seen in Fig. 2, where the CDF of variations against a Gaussian CDF 
is depicted. 



4 Results 

The fluctuation functions F q (s) for timescales ranging from 3 days to iV/5 
are calculated, where iV is the total length of the time series, and for q 
varying between -10 and 10, with a step of 0.5. Fig. 3 shows the MF-DFA2 
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fluctuations F q (s) for various q's. 

A crossover with great magnitude (like as a phase transition) in fluctuation 
function is seen for negative q values in the range 30 < s < 65. The position 
of crossover doesn't have sensitivity to decreasing or increasing q values. 
The only interest behaviour is the asymptotic behaviour of F q (s) at large 
times s. One can clearly observe that above the crossover region, the F q (s) 
functions are straight lines in the double logarithmic plot, and the slopes 
increase slightly when going from high positive moments towards high neg- 
ative moments (from the top to the bottom in Fig. 3). 
For the sake of better studying the large fluctuations, randomized data ( 
both of reshuffled and surrogate data) have been used. Fig. 4 indicates 
that, the magnitude of change in crossover for reshuffled data is very large 
relative to the surrogate data. In fact, one can say that such an effect origi- 
nates mainly from temporal correlations. Moreover, the position of crossover 
is intended to left (about s ~ 4) because of randomizing . 
Monofractal time series are associated with a linear plot r(g), while multi- 
fractal ones possess the spectra nonlinear in q. The highest nonlinearity of 
the spectrums, the strongest multifractality in time series. Calculations in- 
dicate that the time series of exchange rate variations can be of multifractal 
nature. In order to visualize the scaling character of the data, in Fig. 5, the 
corresponding multifractal spectra is shown. Fig. 5 shows three examples of 
r(q) for the original (solid), surrogate (dotted) and reshuffled (dashed) data. 
The nonlinearity of r(q)'s is much weaker for the modified time series than 
for the original ones. Additionally, surrogate data show less nonlinearity 
based on Fig. 5 and therefore, their contribution to multiscaling relative to 
reshuffled data is less. 

The h{q) spectra has been fitted in the range — 10 < q < 10 for original, 
reshuffled and surrogate series by Eq. (10). Representative example for 
original series is shown in Fig. 6. The dotted line in Fig. 6 is obtained by 
best fits of h(q) by Eq. (10). The respective parameters a and b for original, 
reshuffled and surrogate series are listed in Table 2. It is notable that in 
each single case, the q dependence of h(q) for positive and negtive values of 
q can be characterized very well by the two parameters, and all fits remain 
within the error bars of the h(q) values. 

Table 2: Multifractality finger prints (parameters a and b) and strength for 
original, reshuffled and surrogate data. 



Time series 


a 


b 


Aa 


Original data 


0.03 


1.07 


3.54 


Reshuffled data 


0.51 


0.93 


0.60 


Surrogate data 


0.69 


0.8 


0.15 
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It is seen that the strength of multifractality in rial-dollar exchange rate 
variations is very powerful. Moreover, multifractalty strength in randomized 
data decreases specially in surrogate data based on values in Table 2. 
In order to visualizing and better understanding the strength of multifrac- 
tality for original, reshuffled and surrogate data, the singularity spectrum of 
series are shown in Fig. 7. Both the widths of the f(a) spectra in each ran- 
domized data are much smaller than for the original one. This behaviour of 
the reshuffled time series confirms that the persistent autocorrelations play 
an important role in multiscaling of the price variations. But, The spectra 
for the surrogates are typically much narrower than for the reshuffled data 
which can be interpreted as an evidence of the influence of extremely large 
non-Gaussian events on the fractal properties of the time series. 

5 Conclusions 

The multifractal properties of the Iranian rial-US dollar exchange rate log- 
arithmic variations has been studied in this paper through multifractal de- 
trended fluctuation analysis. It is shown that the time series for exchange 
rate variations exhibit the characteristics that can be interpreted in terms 
of multifractality. Its degree expressed by e.g. the widths of the singularity 
spectra f(a) indicate a strong multifractality. Moreover, although the most 
multifractality of the exchange rate variations data is due to different long- 
range correlations for small and large fluctuations, the shape of the proba- 
bility distribution function also contributes to the multifractal behaviour of 
the time series. 



8 



References 

[1] B.B. Mandelbrot, The Fractal Geometry of Nature, Freeman, San Fran- 
cisco, 1983. 

[2] T. Vicsek, Frcatal Growth Phenomena, World Scientific, Singapore, 
1998. 

[3] A-L. Barabasi and H.E. Stanley, Fractional Concepts in Surface 
Growth, Cambridge University Press, New York, 1995. 

[4] T.C. Halsey, M.H. Jensen, L.P. Kanadanoff, I. Procaccia and B.I. 
Shraiman, Phys. Rev. B33, 1141, 1986. 

[5] G. Paladin and A. Vulpiani, Phys. Rep. 156, 147, 1987. 

[6] F. Family and T. Vicsek, J. Phys A18, L75, 1985. 

[7] R. Jullien and R. Botet, Phys. Rev. Lett. 54, 2055, 1985. 

[8] P. Freche, D. Jullien and H.E. Stanley, J. Phys. A18, LI 163, 1985. 

[9] P. Meakin, R. Jullien and R. Botet, Europhys. Lett. 1, 609, 1985. 

[10] J. Lee, Phys. Rev. Lett. 61, 2945, 1988. 

[11] T. Tel, Z. Naturforsch 43a, 1154, 1988. 

[12] J.D. Farmer, Z. Naturforsch 37a, 1304, 1982. 

[13] R. Benzi, G. Paladin, G. Parisi and A. Vulpiani, J. Phys. A17, 3521, 
1984. 

[14] M. Pasquini and M. Serva, Economics Lett. 65, 275, 1999. 

[15] K. Ivanova and M. Ausloos, Physica A 265, 279, 1999. 

[16] A. Bershadskii, Physica A 317, 591, 2003. 

[17] T. Di Matteo, T. Aste and M.M Dacorogna, |cond-mat/0403681| , 2004. 

[18] A. Fisher, L. Calvet and B. Mandelbrot, Multifractality of 
Deutschemark/US Dollar Exchange Rates, Cowles Foundation Discus- 
sion Paper 1166, 1977. 

[19] N. Vandewalle and M. Ausloos, Eur. Phys. J. B 4, 257, 1998. 

[20] A. Bershadskii, Eur. Phys. J. B 11, 361, 1999. 

[21] K. Matia, Y. Ashkenazi and H.E. Stanley, Europhys. Lett. 61, 422, 
2003. 



9 



[22] P. Oswiecimka, J. Kwapien and S. Drozdz, Physica A, 2004. 

[23] H. Takayasu, M.P. Okazaki, K Marumo and T. Shimizu, Fractal Prop- 
erties in Economics, Paradaigm of Complexity, World Scientific, ed. 
M.M Novak (2000). 

[24] F. Schmitt, D. Schertzer and S. Lovejoy, Appl. Stochastic Models Data 
Anal. 15, 29-53 (1999). 

[25] C.K. Peng, S.V. Buldyrev, S. Havlin, M. Simmons, H.E. Stanley and 
A.L. Goldberger, Phys. Rev. E 49, 1685, 1994. 

[26] J. Feder, Fractals, Plenum, New- York, 170, 1988. 

[27] A. Barabasi, T. Vicsek, Phys. Rev. A, 44, 2730, 1991. 

[28] J.W. Kantelhardt, S.A. Zschiegner, E.K. Bunde, S. Havlin, A. Bunde 
and H.E. Stanley, Physica A, 316, 87, 2002. 

[29] J. Theiler, S. Eubank, A. Longtin, J.D. Farmer, Physica D 58 (1992) 
77-94. 

[30] P.F. Panter, Modulatiom, Noise and Spectral Analysis Applied to In- 
formation Transmission, New York, New York, 1965. 

[31] E. Koscielny-Bunde, J.W. Kantelhardt, P. Braun, A. Bunde, S. Havlin, 
Water Resour. Res., |physics/0305078| (2003). 



10 



10000 




500 1000 1500 2000 2500 3000 3500 4000 4500 

t (in days) 



Figure 1: Daily closure rial-dollar exchange rates history (1989-2003). 
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Figure 2: Cumulative distribution function of rial-dollar exchange rate vari- 
ations against a Gaussian cumulative distribution 
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Figure 3: The multifractal fluctuation function F q (s) obtained from multi- 
fractal DFA2 for variations of rial-dollar exchange rates in the period 1989 
to 2003. 
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Figure 4: The multifractal fluctuation function F q (s) obtained from multi- 
fractal DFA2 for randomized (reshuffled and surrogate) variations of rial- 
dollar exchange rates in the period 1989 to 2003. 




Figure 5: Comparison of the miltifractal spectra r(q) of the original and 
randomized exchange rate variations: original (solid), surrogate (dotted) 
and reshuffled (dashed) time series. 
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Figure 6: The generalized Hurst exponents h(q) for the rial-dollar exchange 
rate variations in period 1989 to 2003. The fitted curve has been shown by 
dotted line. 




Figure 7: Comparison of the singularity spectra for original and randomized 
data: original (solid), reshuffled (dotted) and surrogate (dashed) time series. 
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